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1. Introduction 

The investigation of strangeness production in relativistic heavy ion collisions has been 
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proven to be a powerful tool for the study of highly excited nuclear matter, both in terms 
of the reaction dynamics and in terms of its hadrochemistry [HI21IH1IHE1IE10IH1 l^j ITU]. 
Furthermore, strangeness has been suggested as a signature for the creation of a Quark- 
Gluon-Plasma (QGP) PQElEj: in the final state of a heavy- ion collision, subsequent 
to the formation and decay of a QGP, strangeness has been predicted to be enhanced 
relative to the strangeness yield in elementary hadron+hadron collisions. It was actually 
this early work on strangeness as a QGP signature which was instrumental in creating 
the SQM conference series. However, as witnessed by the work presented at this 
conference and its predecessors, strangeness production in heavy-ion collisions has 
5h ' developed into an impressively versatile tool for the characterization of QCD matter. 



In this write-up I will discuss a selection of key contributions at SQM 2006 which I 
consider to have a large impact to the current status of the field of strangeness physics 
or which may have the potential to significantly advance strangeness - or in general 
flavor physics - in the near future. 

I would like to point out that this write-up is not a comprehensive summary covering 
all theory topics discussed at SQM 2006 - for example, there have been excellent mini- 
reviews on strangeness in the nuclear astrophysics context by J. Schaffner-Bielich jTJ 
and S. Reddy jTJ] as well as on the extraction of the nuclear equation of state via near- 
threshold kaon production by J. Aichelin and the extension of the periodic table 
into the strangeness and anti-matter sector by W. Greiner ^3] - the coverage of any of 
these would be beyond the scope of this write-up. 
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2. State of the Art: Lessons and Challenges 

2.1. Hadron Yields and Ratios 

The impressive success of statistical models in describing the (strange) hadron 
abundances and ratios at the CERN/SPS and RHIC [U El EE HHj • and the extracted 
7 S values close to 1 have led to the common conclusion that chemical freeze-out in 
ultra-relativistic heavy-ion reactions occurred very close to - or even at hadrochemical 
equilibrium and that this state most likely has been created by a hadronizing QGP. 
Several studies, exploring the systematics of hadron ratios as function of collision system 
and energy have been presented at this meeting, all corroborating the model's success 
and previous findings [To] . 

However, there remain a number of open questions associated with the success of 
the statistical model, many of which transcend the scope of the model and need to 
be addressed either through direct analysis of data beyond yields and spectra or via 
dynamical non-equilibrium theory approaches based on ab-initio calculations: 

• how did the system achieve chemical equilibrium? Conventional calculations 
based on boost-invariant hydrodynamics with rate-equations for quark production 
[TBJ HZl UH| j pQCD rate-equations [TH] or the Parton Cascade Model [20] all indicate 
that chemical equilibration (and strangeness saturation) cannot be achieved during 
realistic life-times of the deconfined phase. It has been suggested (e.g. in [HUE]) 
that the system would be driven toward and come close to chemical equilibrium 
in the subsequent hadronic phase - a scenario which would help to bridge the gap 
between the calculations indicating insufficient equilibration time in the plasma 
phase and the copious SPS and RHIC data apparently close to chemical equilibrium 
at chemical freeze-out. Recent calculations assuming a hadronizing QGP out 
of chemical equilibrium with subsequent hadronic rescatting have shown that 
rescatting via binary collisions in the hadronic phase is insufficient to drive the 
system toward chemical equilibrium before the expansion of the system leads to 
chemical freeze-out [21J. Currently the most favored approaches for explaining the 
rapid thermalization of the medium are either centered around turbulent color fields 
[2*2*1 123 EH E3 EH] or multi-particle collisions (which may take place either in the 
deconfined or dense, confined phase of the reaction) [23 ESI EH ED]- A smoking 
gun signature for either mechanism remains yet to be established. 

• do the temperature and chemical potential extracted from the statistical model 
fits to final state hadron yields or ratios really reflect the thermodynamic state 
of the system at one particular time during its evolution (i.e. the conditions at 
chemical freeze-out) or are they rather the result of a superposition of different 
states, due to individual hadron species decoupling continuously from the system 
(as would be expected from their different mean free paths)? The latter view is 
supported by a transport model analysis of the time-evolution of the temperature 
and chemical potential in the central cell of a heavy-ion reaction """"Q. However, 
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one would expect a fairly large chi-squared for the single source statistical model 
fit in case of a continuous decoupling scenario, which is not observed in the data. 
It should also be noted that the use of 47r integrated yields at SPS energies creates 
ambiguities in the fit, since many hadron ratios have been shown to exhibit a strong 
rapidity dependence at SPS (and lower) beam energies [H21IHHI- True progress on 
this issue can only be achieved if a model-independent method is found to determine 
the chemical decoupling time of individual hadron species. 

• taking the extracted temperature and chemical potential values for SPS and RHIC 
at face value, the system chemically decouples very close to Tc - at a temperature 
at which the properties of hadrons (e.g. their mass and/or width) could still be 
substantially modified by the temperature and density of the medium (see e.g. the 
spectral function of the p presented at this meeting by A. Foerster [HI])- Is the use 
of hadronic vacuum masses and widths in the statistical model the right approach 
or should their temperature- and density dependent medium modifications be taken 
into account in the model fits (25, 36J ? 

Overall, to this day the physics mechanisms and driving forces behind the impressive 
success of the statistical model are not well understood - finding an answer to the 
question why the statistical model performs so well is a challenge the theory community 
must address with renewed vigor in the near future. 

2.2. Excitation Functions 

Over the past decade the experimental programs at the AGS, SPS and RHIC have 
resulted in a wealth of data which can be compiled into excitation functions of yields, 
ratios, flow coefficients, system sizes etc.. However, the hope of finding some kind of 
sharp discontinuity as a signal of a phase transition has not come to fruition. This may 
be partially due to many QGP signatures predicted to exhibit a sharp discontinuity 
or local extrememum relying heavily on the deconfinement phase-transition being of 
first order with a long-lived mixed phase - an assumption which is at odds with recent 
findings from Lattice QCD, predicting the phase transition to be a continuous crossover 
in the RHIC and SPS domain. Even in the case of a 1st order phase-transition the 
expectation of sharp features in excitation functions may have been somewhat naive: 
sharp discontinuities are predicted as function of temperature, which does not have a 
one-to-one relation to beam energy. Realistic temperature and density profiles as well 
as corona effects (see e.g. the contribution by K. Werner [HI!) w °uld lead to a smearing 
out of the features associated with deconfinement and thus would make for a smooth 
variation in the excitation function. 

Spurred by recent lattice-gauge theory calculations of QCD at finite baryon density 
[HE] there has been renewed hope that an excitation function in incident-beam energy 
might yield interesting results if one can create a QCD medium close to the tri-critical 
point. It has been found that the net-quark susceptibility near the tri-critcial point 
diverges, which would be experimentally accessible via charge fluctuation measurements 
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2.3. Transport Theory 

Relativistic Fluid Dynamics (RFD, see e.g. |HJ H21 EHj ) i s ideally suited for the high- 
density phase of heavy-ion reactions at RHIC, but breaks down in the later, dilute, stages 
of the reaction when the mean free paths of the hadrons become large and flavor degrees 
of freedom are important. The biggest advantage of RFD is that it directly incorporates 
an equation of state as input and thus is so far the only dynamical model in which a 
phase transition can explicitly be incorporated. In the ideal fluid approximation (i.e. 
neglecting off-equilibrium effects) and once the initial conditions for the calculation have 
been fixed, the EoS is the only input to the equations of motion and relates directly 
to properties of the matter under consideration. Ideally, either the initial conditions 
or the EoS should be determined beforehand by an ab-initio calculation (e.g. for the 
EoS via a lattice-gauge calculation), in which case a fit to the data would allow for the 
determination of the remaining quantity. 

RFD has been extremely successful in describing single particle spectra and 
collective flow effects at RHIC O 132]; even though no hydrodynamical model 
implementation has so far attempted to address the entire array of available data in a 
single consistent calculation (a forthcoming publication aims to remedy this situation 
|4"%j). The shape of the spectra as well as the transverse momentum dependence of the 
elliptic flow for minimum bias data are generally reproduced nicely. However, more 
specific centrality bins pose a problem for the elliptic flow calculations. One should 
also bear in mind that hydrodynamical calculations which assume a standard chemical 
equilibrium hadron gas equation of state below T c (which implies simultaneous chemical 
and thermal freeze-out) are unable to fit the measured particle yields simultaneously 
together with the spectra, since statistical models show chemical freeze-out to occur 
around T=170 MeV, whereas the shape of the spectra requires a hydrodynamic evolution 
to T=110 MeV. One method to deal with the separation of chemical and thermal freeze- 
out is the partial chemical equilibrium model (PCE) [1§J E3 OH]- Below a chemical 
freeze-out temperature T c h one introduces a chemical potential for each hadron whose 
yield is supposed to be frozen out at T c h. The PCE approach can account for the proper 
normalization of the spectra, however, it fails to reproduce the transverse momentum 
and mass dependence of the elliptic flow. More importantly for the strangeness sector, 
ideal RFD calculations lack the capability of dealing with the flavor dependence of 
hadronic cross sections and a possible flavor-dependent sequential freeze-out: recent 
experimental results suggest that at thermal freeze-out multistrange baryons exhibit less 
transverse flow and a higher temperature closer to the chemical freeze-out temperature 
compared to non- or single-strange baryons IHH] . This behavior can be understood 
in terms of the flavor dependence of the hadronic cross section, which decreases with 
increasing strangeness content of the hadron. The reduced cross section of multi-strange 
baryons leads to a decoupling from the hadronic medium at an earlier stage of the 
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reaction, allowing them to provide information on the properties of the hadronizing 
QGP less distorted by hadronic final state interactions jSH 153]. 

The reach of RFD can be extended and the problem having to terminate the 
calculation at a single flavor-independent and fixed freeze-out temperature can be 
overcome by combining the RFD calculation with a microscopic hadronic cascade model 
- this kind of hybrid approach (dubbed hydro plus micro) was pioneered in [SH] and 
has subsequently been adopted by other groups [SHI ED EH EE]- Its key advantages 
are that the freeze-out occurs naturally as a result of the microscopic evolution and 
that flavor degrees of freedom are treated explicitly through the hadronic cross sections 
of the microscopic transport. Due to the Boltzmann equation being the basis of the 
microscopic calculation in the hadronic phase, viscous corrections for the hadronic phase 
are by default included in the approach - the full treatment of viscosity in the deconfined 
phase in a 3D hydrodynamic calculation remains a challenge for the future. 

Such hybrid macro/micro transport calculations are to date the most successful 
approaches for describing the soft physics at RHIC [H2]- First implementations of 
these hybrid approaches were restricted to 1 + 1 jSH] and 2+1 jJH] dimensions in the 
hydrodynamic component of the model. However, recently new state of the art fully 
three-dimensional ideal hydrodynamics models have become available and are now being 
incorporated into the hydro+micro framework [17| IHUj. loT} I48j . 

The latest model implementation has been shown at this conference [IB] to be 
able to excellently reproduce the measured Pt spectra and (Pt) values of hyperons and 
multi-strange baryons at RHIC (see left frame of figure 1). A collision number analysis 
for these hadron species based on this calculation [IS] confirms the early and sequential 
freeze-out findings described in 155]. 

Marcus Bleicher discussed elliptic flow at RHIC in the framework of the microscopic 
hadron/string UrQMD and RQMD transport models jHSUHHl- Both models reach only 
60% of the absolute magnitude of the measured t>2, due to the string degrees of freedom 
being used to describe the initial phase in these models generating insufficient pressure. 
Interestingly, the model calculations exhibit features reminiscent of constituent quark 
number scaling in the intermediate px range between 2 and 6 GeV/c. This is due to the 
use of the Additive Quark Model for most hadron-hadron cross section implemented in 
the model, thus having the interaction scale with the number of constituent quarks. At 
closer inspection, however, one finds a flavor ordering in the elliptic flow generated by 
the microscopic transport models: Vz(N) > V2((Y) > t>2(5) > t>2(fi), due to the strange 
quark being assigned a smaller interaction cross section than the u and d quarks. This 
structure is not observed in the data, which rather indicate the build-up of elliptic flow 
independent of the quark flavor prior to hadronization. 

2.4- Hadronization via Parton Recombination 

The detailed understanding of hadronization plays a crucial role for isolating signatures 
sensitive to the QGP evolution and properties of the system from those which are 
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Figure 1. Left: Pt spectra for A, 5, fl in central collisions calculated in the hybrid 
3D-hydro+UrQMD approach compared to STAR data. Right: elliptic flow of (multi- 
strange hadrons measured by the STAR collaboration (SHEII- Note the <j> meson 
following the kaon elliptic flow. 



dominated by the later reaction stages. One of the theoretical milestones of the first 
several years of the RHIC program was the development of the recombination plus 
fragmentation model as the standard model of hadronization for matter at intermediate 
and high transverse momenta: at the center of the recombination + fragmentation model 
is the realization that hadron production at momenta of a few GeV/c in an environment 
with a high density of partons occurs by recombination, rather than fragmentation, of 
partons. It is found that recombination always dominates over fragmentation for an 
exponentially falling parton spectrum, but that fragmentation wins out eventually, when 
the spectrum takes the form of a power law. 

Among the RHIC discoveries that prompted the development of the recombination 
models of hadronization, is that the amount of suppression of hadron production at 
intermediate px compared to the scaled proton-proton baseline seems to depend on the 
hadron species. In fact, in the production of protons and antiprotons between 2 and 
4 GeV/c the suppression seems to be completely absent. Generally, pions and kaons 
appear to suffer from a strong energy loss while baryons and antibaryons do not. Two 
stunning experimental facts exemplify this |63[ IMj . First, the ratio of protons over 
positively charged pions is equal or above one for px > 1.5GeV/c and is approximately 
constant up to 4 GeV/c. Second, the nuclear suppression factor Raa below 4 GeV/c is 
close to one for protons and lambdas jH7J, while it is about 0.3 for pions. 

The recombination approach jHHl EM EH HI] has been able to account for the above 
baryon/meson differences. Additionally, it was observed that the elliptic flow pattern 
of different hadron species can be explained by a simple recombination mechanism 
[721 H3 EH US]- The anisotropies ,t> 2 , for the different hadrons are compatible with 
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a universal function ^(pt) i n the parton phase, related to the hadronic flow by factors 
of two and three depending on the number of valence quarks [Tol I77j . 

One of the perceived weaknesses of the recombination+fragmentation approach was 
that its development was triggered by the experimental observation of the meson/baryon 
anomalies in the RHIC data and therefore its success in explaining these features 
occurred after the fact. However, the recombination approach was able to predict the 
the elliptic flow t> 2 of multi-strange hadrons such as the 0, 5 and Q as a function 
of transverse momentum |78j : the measurement of V2 for the and Q allows for the 
unambiguous distinction between parton recombination and statistical hadro-chemistry 
to be the dominant process in hadronization at intermediate transverse momenta, since 
e.g. the meson has approximately the mass of a nucleon, but the valence-quark 
content of a meson. In a hydrodynamic picture with the hadron mass as the guiding 
scale the would follow the systematics of the nucleon whereas in the recombination 
picture it would follow the behavior of the pions and kaons. Data on the elliptic flow 
of the meson and Q baryon presented at this conference, e.g. in the talk by S.L. 
Blyth - see the right frame of figure 1 - clearly exhibit constituent quark scaling and 
therefore impressively confirm the physics of the recombination model and demonstrate 
its predictive power. 

3. New Directions 

3.1. Quarks in the Color Glass 

Heavy-ion collisions at ultra-relativistic energies can be described by the collision of two 
coherent sheets of high energy density gluonic fields (commonly referred to as Color 
Glass Condensate) [EH EDI EI]- Since the physical density of gluons becomes large, 
their typical separation is small, implying a small value for a s . Furthermore, these 
highly coherent gluons saturate the phase space up to the maximal occupation number 
~ l/a s . Due to the weak coupling, it is possible to describe this system from first 
principles in QCD. 

The Color Glass Condensate (CGC) has been suggested to describe the initial state 
of gold nuclei in RHIC collisions. While there is a broad consensus that at sufficiently 
high beam energies the saturation physics of the CGC should dominate the initial state, 
it is still a matter of debate whether these conditions are actually fulfilled at RHIC. 
Nonetheless, the further development of the CGC picture for heavy-ion collisions is of 
great theoretical importance. Original work on the CGC focused on the description 
of the initial state as gluonic field. Due to the large gluon densities at low Bjorken-x, 
these are thought to dominate the dynamics of the initial state. However, the lack of 
treatment and consideration of quark and anti-quark production makes it difficult to 
connect the CGC to experimental data, e.g. in the strangeness sector. 

T. Lappi presented a numerical integration of the Dirac equation in order to 
calculate the number of quark- ant iquark pairs initially produced in the classical color 
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fields of colliding ultra-relativistic nuclei (the resulting state of high energy and density 
matter being termed the glasmd) [H2]- While the number of qq pairs is parametrically 
suppressed in the coupling constant, he found that in the CGC their production rate 
is comparable to the thermal ratio of gluons/pairs = 9iVy/32. After isotropization one 
thus would end up with a quark- gluon plasma in chemical equilibrium. This finding 
is of great significance since it provides a link between the CGC initial state and an 
abundance of flavor-centric data indicating the formation of a QGP in thermal and 
chemical equilibrium. 

3.2. Conserved Charge Correlations 

The elementary electric charge, strangeness and baryon number carried by the QGP 
degrees of freedom differs significantly from those in a gas of hadrons. This observation 
does not only apply to the elementary charges themselves, but also to the size of the 
average fluctuations of net baryon number, strangeness and electric charge in a finite 
volume. While hadronization and confinement prohibit us from directly observing the 
fractional electric charge and baryon number and strangeness of the QGP degrees 
of freedom, these event by event fluctuations may under certain conditions survive 
hadronization and subsequent hadronic rescattering and thus serve as indicators of the 
existence of a QGP (SHUSH- 

Surprisingly, all experimental analysis of net charge fluctuations for SPS and RHIC 
data agree with the hadron gas prediction, giving no indication at all about a possible 
deconfined phase. Since many other measurements are compatible with the assumption 
of deconfinement, it is tempting to speculate that the dynamics of hadronization and/or 
hadronic final state interactions strongly affect the charge fluctuation observable and 
mask the fluctuations generated in the deconfined phase. 

A more robust approach of exploiting the difference in conserved elementary 
charges between QGP and HG degrees of freedom lies in the area of correlations: in 
particular, the correlation between baryon number B and strangeness S may elucidate 
the microscopic structure of the QGP and the nature of its degrees of freedom jSnHSHllSZI- 
Strangeness is carried exclusively by s and s quarks which carry baryon number as well: 
B s = —^S, whereas in a hadron gas the relation between B and 5* is much less intimate. 
This Cbs correlation can be directly extracted from Lattice QCD calculations and is 
experimentally accessible via event-by-event fluctuations: 



For a system of hadrons Cbs can be formulated in terms of multiplicity variances 
a k = (n\) — (n k ) 2 « (n k ) for the respective hadron species k: 



Cbs = ~ 3 
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Figure 2. Left: Cbs and Cqs in a truncated hadron gas at \jlb = MS = MeV 
compared to lattice calculations at fi = 0. The two bands denoting the hadron 
gas calculation and reflect the uncertainty in the actual value of the phase transition 
temperature Tc, which is assumed to lie in the range of 170±f MeV (figure taken from 
86 ). Right: elliptic flow t>2 of non-photonic electrons in central Au+Au collisions at 
RHIC. PHENIX data are compared to theory predictions using Langevin simulations 
with elastic c and b quark interactions in an expanding QGP (figure taken from 



For an ideal gas of hadrons at T = 170 MeV and zero chemical potential hb — one 
finds C BS = 0.66. 

The Cbs correlator can also be expressed in terms of susceptibilities, Cbs — 
—3xbs I Xss which are second derivatives of the free energy with respect to the chemical 
potential and can be directly be calculated from lattice QCD: 

1 d 2 F 1 d 2 F 

XBS — -77 o a J ~X,SS 77 7^ q~ I 3 ) 

V arsons V <9/4 

Using values of lattice susceptibilities extracted at T = 1.5 Tc in [HE] on finds Cbs ~ 1 
indicating that the degrees of freedom of a QGP carry baryon number and strangeness 
of individual quarks and suggesting that the quark flavors are uncorrelated as in an ideal 
QGP (see left frame of figure 2). Note that the presence of pure gluon clusters cannot 
be ruled out by this diagnostic. This finding is of great importance, since the evaluation 
of Cbs f° r different models of the QGP structure allows for the direct verification or 
falsification of these models when compared to data and lattice QCD. 

3.3. Anomalous Viscosity 

As stated earlier, measurements of the anisotropic collective flow of hadrons emitted in 
non-central collisions of heavy nuclei at the Relativistic Heavy Ion Collider (RHIC) are 
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in remarkably good agreement with the predictions of ideal relativistic fluid dynamics. 
In order to describe the data, calculations need to assume that the matter formed in 
the nuclear collision reaches thermal equilibrium within a time T\ < 1 fm/c [89J and 
then expands with a very small shear viscosity r\ <C s, where s is the entropy density 
PU] . The comparison between data and calculations indicates that the viscosity of the 
matter cannot be much larger than the postulated lower bound 77 min = s/47r [HI], which 
is reached in certain strongly coupled supersymmetric gauge theories fiTH. 

This result is nontrivial because the shear viscosity of a weakly coupled, 
perturbative quark-gluon plasma is not small. In fact, the perturbative result for the 
shear viscosity, in leading logarithmic approximation, is [OB] 

g^lng 1 

where df ~ 0(100) is a numerically determined constant that weakly depends on the 
number of quark flavors rif. The result (J3J), as well as the finding that numerical 
solutions of the Boltzmann equation exhibit fluid dynamical behavior only when the 
cross section between gluons is artificially increased by a factor ten or more [01], have 
invited speculations that the matter produced at RHIC is a strongly coupled quark- 
gluon plasma (sQGP). The possible microscopic structure of such a state is not well 
understood at present E3 EH] • 

However, as Berndt Miiller pointed out in his talk [98J, there exists an alternative 
mechanism that may be responsible for a small viscosity of a weakly coupled, but 
expanding quark-gluon plasma. This mechanism is based on the theory of particle 
transport in turbulent plasmas |HZ1 EH]- Such plasmas are characterized by strongly 
excited random field modes in certain regimes of instability, which coherently scatter the 
charged particles and thus reduce the rate of momentum transport. The scattering by 
turbulent fields in electromagnetic plasmas is known to greatly increase the energy loss 
of charged particles |100| and reduce the heat conductivity |101| 1102] and the viscosity 
103, 1104] of the plasma. Following Abe and Niu |104j . the contribution from turbulent 
fields to transport coefficients was called "anomalous" . 

The sufficient condition for the spontaneous formation of turbulent, partially 
coherent fields is the presence of instabilities in the gauge field due to the presence of 
charged particles. This condition is met in electromagnetic plasmas with an anisotropic 
momentum distribution of the charged particles |105j . and it is known to be satisfied 
in quark-gluon plasmas with an anisotropic momentum distribution of thermal partons 

[221 ESI EU - 

The additional contribution to the viscosity, r] B , induced by the turbulent fields 
decreases with increasing strength of the fields. Since the amplitude of the turbulent 
fields grows with the magnitude of the momentum anisotropy, a large anisotropy will 
lead to a small value of rjs- Because the relaxation rates due to different processes are 
additive, the total viscosity is given by 

r j - l = r lB 1 + 7 lc 1 . (5) 
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This equation implies that T]b dominates the total shear viscosity, if it is smaller than 
7]c- In that limit, the anomalous mechanism exhibits a stable equilibrium in which 
the viscosity regulates itself: The anisotropy grows with 77, but an increased anisotropy 
tends to suppress r\ B . The result is that in the weak coupling limit, the anomalous 
viscosity is much smaller than the viscosity due to collisions among thermal partons. By 
reducing the shear viscosity of a weakly coupled, but expanding quark-gluon plasma, this 
mechanism could possibly explain the observations of the RHIC experiments without 
the assumption of a strongly coupled plasma state. 

4. From SQM to F(lavor)QM 

4-1. Heavy-Quark Production, Diffusion and Energy Loss 

The calculation of heavy quark production in the framework of pQCD is fairly well 
established - recent improvements having been reported by I. Vitev at this meeting 
|107j (see also |106| and references therein). Among the novelties presented was the 
work by K. Tuchin, who discussed heavy quark production in High Parton Density 
QCD in a quasi-classical approximation, including low- 2 quantum evolution, as well as 
heavy-quark production based on the effect of pair production in external fields |l()8j . 

For moderate transverse momenta (jp? < a few GeV/c) the energy loss of heavy 
quarks is thought to be dominated by inelastic collisions with medium constituents, 
rather than gluon radiation |1U9|.I110] . because the heavy quarks are not ultra-relativistic 
and gluon radiation is suppressed by the so-called dead-cone effect |lllj . These findings 
were highlighted in a comprehensive analysis of the non-photonic electron nuclear 
suppression factors measured at RHIC presented by M. Djordjevic |112t I113j . The 
analysis clearly indicates the importance of including charm and bottom quarks in the 
calculation and taking their radiative as well as collisional energy loss into account. 

The importance of the collisional energy loss contribution and the possibility of 
heavy quarks actually thermalizing in the medium provides an opportunity to utilize 
them as probes for the transport coefficients of the QCD medium. Since collisional 
energy loss occurs in many small steps, the motion of a heavy quark can thus be described 
by a Fokker-Planck or, equivalently, Langevin equation. Several studies based on such 
an approach have recently been done |114[ II 1^ lllfij . The most detailed and extensive 
one of these was performed by Moore and Teaney |114j . who derived an expression 
for the diffusion coefficient D in the framework of hard-thermal loop (HTL) improved 
perturbation theory and discussed the limitations of the Fokker-Planck approach. These 
authors also studied the resulting phenomenology of heavy quark transport in dense 
matter created by a heavy ion collision, using a boost invariant hydrodynamical model 
with an ideal equation of state. At SQM 2006, heavy quark diffusion calculations were 
discussed by P.B. Gossiaux |116| I117j as well as by R. Rapp |115| I118j . Whereas the 
philosophy of I117j is to determine the transport coefficients from a comparison to 
data, the approach of |115[ I118j is based on introducing a resonant charm - light quark 
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interaction and then calculating the drag- and diffusion coefficients for the Langevin 
evolution on the basis of that interaction. A comparison to data would thus yield 
information on the microscopic interaction between charm and the light quark species 
(see right frame of figure 2). 

4-2. Charmonium Spectral Functions 

In 1986 Matsui and Satz proposed |119j that the suppression of heavy quarkonia-mesons 
could provide one of the signatures for deconfinement in QCD at high temperatures. The 
idea was based on an analogy with the well known Mott transition in condensed matter 
systems. At high densities, Debye screening in a quark-gluon plasma reduces the range 
of the attractive force between heavy quarks and antiquarks, and above some critical 
density screening prevents the formation of bound states. The larger bound states were 
expected to dissolve before the smaller ones as the temperature of the system increases. 
The ip' and Xc states were thus expected to become unbound just above T c , while the 
smaller ip state would only dissolve above ~ 1.2T C . 

However, QCD lattice gauge-theory calculations of charmonium correlators in 
recent years have necessitated a revision in our understanding of the dissociation of 
charmonium states and what it implies for the properties of the surrounding medium: 
as discussed by P. Petreczky at this meeting, the IS charmonia states {J/ip, ip' and 
r/c) survive to unexpectedly high temperatures above 1.5 T c and only the IP states 
(i.e. the \c) dissolves around 1.2 Tc |12()| 1121] . It is therefore questionable whether 
the observed charmonium suppression in Pb+Pb collisions at the SPS |122j is truly 
a smoking gun signature for deconfinement or rather the result hadronic dissociation 
(see e.g. P.23J ) . Moving on to RHIC energies, the situation will become even more 
complicated: achievable temperatures may be above the threshold for charmonium 
dissociation (which is unlikely in the SPS case), however, the additional suppression may 
be compensated through novel production mechanisms. The multiplicity of produced 
charm and anti-charm quarks per event at RHIC is sufficiently large that charmonium- 
regeneration via parton recombination may occur |124[ I125j , giving rise to a fairly flat 
behavior of the observed charmonium yield as a function of beam energy. 

Another exciting prospect is the measurement of charmonium elliptic flow, 
first shown in form of preliminary data by NA60 here at this conference [34J. 
The measurement of significant charmonium elliptic flow would indicate charm 
thermalization and charmonium regeneration via parton recombination. In addition it 
has been shown |118j . that this observable exhibits a strong sensitivity to the in-medium 
interaction of charm quarks. 

5. Outlook 

The field of strangeness in heavy-ion collisions - and with it this conference series - has 
developed tremendously over the past 20 years. While strangeness as a QGP signature 
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was instrumental in creating the SQM conference series, strangeness production in 
heavy-ion collisions has become an impressively versatile tool for the characterization of 
all aspects of confined and deconfined QCD matter, as witnessed by the work presented 
at this conference. The future of the field will most likely bring the generalization of 
the concepts and lessons learned from strangeness into the heavy-quark sector at higher 
incident beam energies. The application of flavor-dominated physics phenomena as 
probes of the the hot and dense QCD medium will be an exciting field for many years 
to come. 
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